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ABSTRACT 

We present the discovery of pulsational variations in the cool magnetic Ap star HD 1 15226 - the first high- amplitude rapidly oscillating 
Ap star discovered with time-series spectroscopy. Using high-resolution spectra obtained with the HARPS instrument at the ESO 3.6- 
m telescope, we detect radial velocity variations with a period of 10.86 min in Prm, Ndm, Dym lines and in the narrow cores 
of hydrogen lines. Pulsational amplitudes exceed 1 kms~' in individual lines of Ndm. The presence of running waves in the stellar 
atmosphere is inferred from a phase shift between the radial velocity maxima of rare-earth and hydrogen lines. Our abundance analysis 
demonstrates that HD 1 15226 exhibits typical roAp spectroscopic signature, notably ionization anomaly of Pr, Nd and Dy. We discuss 
the discovery of pulsations in HD 115226 in the context of recent spectroscopic studies of roAp stars and point to the existence of 
correlation between spectroscopic pulsational amplitude and the stellar rotation rate. 
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1. Introduction 

Rapidly oscillating Ap (roAp) stars are cool magnetic chem- 
ically peculiar stars, pulsating in high-overtone non-radial p- 
modes with periods around 10 min. Excitation of oscillations 
in these stars is closely related to the presence of strong global 
magnetic field, as the pulsational axis is seen to be aligned with 
the axis of oblique magnetic field and pulsations are enhanced at 
the magnetic poles (Kochukhov 2004). In addition to remarkable 
pulsation behav iour, which opens interes ting prospects for as- 
teroseismology d Vauclair & T heado 2004), roAp stars are char- 
acter ized by extremely anoma lous surface chemical composi- 
tion (Ryabchi kova et all 120041) . as well as strongly inhomoge- 
neous horizontal and vertical distribution of chem ical elements 
dKochukhov et al.ll2004tlRvabchikova et al.ll2002l) . 

Most of the previous work on identification and frequency 
analysis of the three dozen currently known roAp stars was 
performed using time-s eries photometry at small telescopes 
dKurtz & Martinezi2000i) . However, inspired by the discovery of 
outstanding pu lsational variations in the lines of rare-earth ele- 
ment s (REEs) dSavanov et al . 1999t lKochukhov & Ry abchikova 
1200 lb . the focus of observational studies of roAp stars is shifting 
towards ti me-resolved spectroscopy. A number of recent stud- 
ies (e.g jMkrtichian et al.ll2003l:lRvabchikova et al.ll2007bl) have 
demonstrated unique diagnostic potential of the precise radial 
velocity and line profile measurements in individual lines. This 
type of observational material allows one to trace propagation 
of pulsation waves through the chemically stratified stellar at- 
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mosphere, eventually obtaining a three-d imensional picture of 
pulsations and chemical inhomogeneities dKochukhovl2007t3) . 

Despite these noteworthy achievements of time-resolved 
spectroscopic studies, relatively little is done in the direction of 
expanding the list of known roAp stars and understanding the 
coexistence of pulsating and apparently constant Ap stars in the 
same region of the H-R diagramm. With the aim to obtain an 
improved picture of the incidence of pulsations in Ap stars, we 
have initiated a survey of cool peculiar magnetic stars with the 
ultra-stable spectrograph HARPS at ESO. Here we report the 
first result of our observations - the discovery of high-amplitude 
oscillations in the magnetic Ap star HD 1 15226. 

2. Observations and data reduction 

We have observed HD 115226 on the night of Ap ril 15, 2007, 
with the HARPS spectrograph dMavor etal .112003b at the ESO 
3.6-m telescope at La Silla. The observations have started at the 
barycentric JD 2454205.63472 and continued for 4.3 hours. In 
total, 102 consecutive 120 s exposures were collected. The dead 
time between stellar observations was 31s. 

The extraction of one-dimensional spectra and barycentric 
velocity correction of the wavelength scale was carried out with 
the HARPS pipeline available at the telescope. The peak signal- 
to-noise ratio of 40-45 was achieved in the individual expo- 
sures of HD 1 15226. The spectra have resolving power A/AA = 
115 000 and cover a wavelength range from 3780 to 6910 A, 
with a 30 A gap near 5320 A. 

The dominant source of radial velocity noise in our data will 
be due to the photon noise rather than to the instrumental pre- 
cision. Therefore, we did use the simultaneous ThAr method, 



avoiding in this way any contamination of the stellar signal. At 
the beginning and at the end of HD 1 15226 observations we took 
a ThAr reference spectrum, and can estimate that the instrumen- 
tal drift within the time series was at most 0.3 6 ms 4 . 

Fo llowing the procedure outlined in iKochukhov et al.l 
(2007b), we post-processed all extracted spectra of HD 115226 
with the aim to achieve an accurate and consistent continuum 
normalization. 

3. Basic properties of HD 11 5226 

Little is known about HD 115226 (HIP 64883, CP-72 1373). 
This southern V = 8.51 chemically peculiar star is clas- 
sified as A 3p Sr in th e Gen er al Catalogue of Ap and 
Am Stars dRenson et all 1 199 lb . iGomez et ail ([998) and 
Koch ukhov & Bagnulol (120061) included this object in their stud- 
ies of the evolutionary state of magnetic Ap stars. Using 
Hipparcos parallax and T e fi = 7640 K estimated from the Geneva 
photometric colors, the authors of the latter paper determined 
L = 7.2 ± 1 .9 Lq. F rom the comparison with standard stellar evo- 
lutionary models (ISchaller et al.ll 1 992b one can infer that these 
parameters correspond to a M — 1 .60 + 0.05 M Q star with a frac- 
tional age between 0.0 and 0.4 of the main sequ ence lifetime. 

The Stromgren photometry of HD 1 15226 (iMartinezll 1993b 
indicates T e g = 8000-8200 K i f normal-star cal i bration s of 
iMoon & Dworetskyl (1 19851) and iNapiwotzki et al.l (1 19931) are 
used. However, comparison of the theoretical calculations with 
the observed profiles of Ha and H/3 favours lower r e ff suggested 
by the Geneva photometric param eters. 

IKochukhov & Bagnulo (2006) determined longitudinal field 
(B z ) = 752 + 48 G from a single circular polarization measure- 
ment of HD 115226 with the FORS1 instrument at VLT. This 
robust magnetic field detection, along with the typical peculiar- 
star appearance of the spectrum of HD 115226 (Fig. [TJ, leaves 
no doubt that this star belongs to the group of cool magnetic Ap 
stars. The value of inferred for HD 1 15226 is also typical of 
the known roAp stars. 

To perform a preliminary analysis of the atmospheric chem- 
istry of HD 115226 and to compare it with the chemical abun- 
dances of well-studied roAp stars, we adopted a model atmo- 
sphere of the roAp star 10 Aql (HD 176232) with the parameters 
T e s = 7650 K, lo g # -4.0, determined in ou r study of Ca stratifi- 
cation in 10 Aql (Ryabc hikova et ai1l2007al) . A rather rapid rota- 
tion of HD 1 15226 does not allow us to study individual lines for 
most elements, therefore we perform ed magnetic spectra l syn- 
thesis with the help of synthmag code (lKochukhovl2007bl) . First 
of all, we tried to estimate v e sin i and the strength of surface 
magnetic field using magnetic intensification. Line profiles of 
almost all elements in the spectrum of HD 115226 have very 
complex structure, indicating a high degree of the horizontal 
abundance inhomogeneity. Best-fit projected rotational veloci- 
ties vary from 18 kms~' for Coi and Yn to 23-25 krns" 1 for 
REEs and Si n lines and up to 28-30 km s _1 for Mg, Ca, Cr and 
Fe lines. This dispersion in broadening is accompanied by dif- 
ferent shifts of the line centroids: -10 kms 1 for Coi, from -5 
to -8 kms 1 in REEs, from -2 to -3 kms 1 in Mg, Ca, Cr, 
Fe, and +7 kms -1 for Si. These discrepancies point to a spotted 
distribution of chemical elements at the stellar surface. 

Our spectrum synthesis was performed for the magnetic field 
strength in the range from 1 to 3 kG. The lowest dispersion in the 
abundances inferred from spectral lines with different Lande fac- 
tors (0-2.5 for Fe lines and 0.6-1.6 for Ndm lines) and different 
Zeeman structure was obtained for the magnetic field strength of 
1 kG. 
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Fig. 1. Comparison of the 6120-6170 A region in the average 
spectrum of H D 1 15226 with the know n rapidly rotating roAp 
star HD 83368 (IKochukhov et al.ll2004 . Prominent lines of Prm 
and Ndm appear in both stars. The spectrum of HD 83368 is 
shifted downwards for display purpose. 

Table 1. Chemical abundances in the atmosphere of HD 1 15226. 
Error estimates are given for species with sufficient number of 
measured lines. In other cases we provide a range of abundance 
that fits individual spectral features. 



Ion 


log(2V e i/tftot) 


Ion 


log(AWAU 


Mg i/ii 


-4.5 


Cen 


-9.5 


Si i/ii 


-4.6 


Pru 


-9.6 


Cai/n 


-5.4 ±0.2 


Prm 


-8.1 ±0.2 


Tin 


-7.0 


Ndn 


-8.8 to -9.0 


Cri/n 


-5.0 


Ndm 


-7.2 ±0.3 


Fe i/ii 


-4.5 to -4.6 


Smu 


-9.2 to -9.5 


Nil 


-6.3 


Eun 


-10.2 


Coi 


-6.0 


Gdn 


-9.3 


Yn 


-8.5 


Tbm 


-8.2 


Ban 


-10.5 


Dyn 


-9.0 


La n 


-9.50 


Dy m 


-7.3 ±0.2 



The result of abundance analysis of HD 115226 is summa- 
rized in TableQ] The inferred abundances are typical of roAp 
stars in the 7500-8000 K effective temperature range. The star 
has nearly solar concentration of Mg, Si, Ti, Fe, overabundance 
of Cr and Co, underabundance of Ba. The well-known REE 
anomaly, where abundance derived from the lines of second ions 
(Prm, Nd in, Dy m) exceeds by 1 .5-2 dex those derived from th e 
lines of singly ionized species (Ryabc hikova et af]|200ll 12004). 
is also present. Interestingly, the vertical Ca abundance distribu- 
tion determined for the atmosphere of 10 Aql fits well the Ca 
line profiles in the spectrum of HD 1 15226. This suggests close 
similarity of the Ca stratifications in the atmospheres of these 
two stars. 

Significant Doppler broadening of the HD 1 15226 spectrum 
indicates that the stellar rotational period is short. Combining 
v e sin/ = 25-30 krns -1 determined here with the information 
on stellar luminosity and temperature, we find P Iot < 3.0-3.5 d. 
However, despite the presence of surface abundance inhomo- 
geneities, HD 1 15226 does no t exhibit a prominent variability 
in the ASAS (|Pojmanskill2002l) photometric time-se ries. At the 
same time, Hipparcos Epoch Photometry of this star (lESAl 19971) 
shows a marginal variation with P — 3.61 d and an amplitude be- 
low 0.01 mag. 



4. Analysis of pulsational variability 

Large rotational broadening in the spectrum of HD 1 15226 sig- 
nificantly reduces the number of lines suitable for precise ra- 
dial velocity measurements. Fortunately, the star shows strong 
lines of Prm and Ndm, which e xhibit pulsational variability 
in almost all known roAp stars dKochukhov & Ryabchikova 
1200 It IRvabchikova et all l2007bl) . Using the VALD database 
(lKupkaetal.il 19991) an d an i mproved Ndm line list published 
by IRvabchikov a et al.l ([2006), we measured radial velocities of 
the absorption features of doubly ionized REEs and the hydro- 
gen line cores. Velocities were determined with the centre-of- 
gravity method. Based on these measurements, the amplitudes 
and phases of variable lines were found by least-squares fitting 
of the cosine function A cos [2n(t/P + <p)]. 

The radial velocity analysis of strongest rare-earth lines read- 
ily shows pulsational signal with the frequency w 1.53 mHz. The 
cores of Ha and H/3, as well as weaker lines of doubly ionized 
Pr, Nd and Dy, also show variation with this frequency. Thus, 
our observations demonstrate that HD 115226 is a new roAp 
star, pulsating with a period typical of this class of oscillating 
stars. However, the radial velocity amplitude in HD 1 15226 sig- 
nificantly exceeds the average level of few hundred ms~' seen 
in other roAp stars. 

In order to determine pulsational characteristics of 
HD 115226 with a better precision, we analysed combined 
radial velocity measurements of H 1 (2 lines), Pr m (5 lines) and 
Ndm (17 lines). Corresponding amplitude spectra are presented 
in Fig. |2] Combined velocity of the two weak Dy m lines shows 
higher noise level (lower panel in Fig.|2|, therefore we disregard 
this REE ion in the period determination. The least-squares fit 
yields P = 10.880±0.031,10.852+0.016andl0.867+0.010min 
for H 1, Pr m and Nd 111, respectively. The weighted mean value 
of the pulsation period is thus 10.864 + 0.008 min (frequency 
1.534 mHz). Fig. [3] shows the average radial velocity curve 
of Ndm together with a sinusoid fit. Analysis of the residuals 
shows no evidence of additional frequencies with amplitudes 
above 110 ms 1 . For Hi and Prm the corresponding noise 
level (amplitude of the highest noise peak in the residuals) 
is 160 ms~'. Further spectroscopic observations at a larger 
telescope are required to confirm monoperiodic character of 
oscillations in HD 1 15226. 

Finally, we have derived least-squares estimate of the ampli- 
tudes and phases of 26 variable lines for the fixed 10.864 min 
period (Table [2]). Integration time of 120 s adopted in our ob- 
servations of HD 115226 is non-negligible compared to the in- 
ferred pulsational period. We found that for the time sampling 
of our data the phase smearing reduces amplitudes by 6.2 %. 
Consequently, all radial velocity amplitudes reported here are 
multiplied by a factor 1.066 to correct for this effect. 

From Table [2] it is evident that oscillation amplitudes show 
considerable scatter even for the lines of the same ion. This can 
be attributed to the dilution of pulsation signal in the lines af- 
fected by unrecognized blends. For this reason our amplitude 
determination possibly gives only a lower limit of the real signal 
for many pulsating lines. Thus, we believe that the radial veloc- 
ity amplitudes exceeding 1 km s _1 , as seen for Nd m 4914, 5802, 
5851, 6327 and 6550 A, are more representative of the pulsa- 
tions in HD 1 15226 than variability in the lines with amplitudes 
around 500 ms~'. 

Time-series analysis reveals that the lines of Pr m and Nd m 
pulsate with approximately the same phase, which differs from 
that of Dy m and H 1. In our convention a larger pulsation phase 
corresponds to an earlier radial velocity maximum. Therefore, 
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Fig. 2. Amplitude spectra for the combined radial velocity of two 
Hi 5 Prm, 17 Ndm, and two Dy m lines. The presence of rapid 
oscillation with the frequency 1.53 mHz is evident. 
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Fig. 3. Average radial velocity of 17 Nd m lines as a function of 
time. Symbols show observations, the solid curve illustrates a 
least-squares cosine fit with P = 10.864 min. 



we find that the Dy m velocity lags by 0. 1 period behind H 1. The 
phase lag reaches 0.2 for Pr m and Nd m. 



Table 2. Amplitudes and phases of individual lines of H i, Pr m, 
Nd m and Dy in determined by the least-squares analysis of ra- 
dial velocity measurements. The phases, determined with re- 
spect to the BJD of the first observation (see Sect. 13, are given 
in units of pulsation period. 
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A number of strong lines belonging to the iron-peak ele- 
ments and Ca were also searched for pulsation signatures. No 
variability exceeding 100 m s~' was detected in these lines. 

5. Discussion 

We show that HD 115226 is a new roAp star with broad lines 
and prominent pulsational variability in the cores of hydrogen 
lines and in doubly ionized Pr, Nd and Dy. The radial velocity 
amplitudes, exceeding 1 kms 1 in some Ndm lines, are rather 
high for a roAp star. On the other hand, atmospheric parameters 
and chemical abundance pattern of HD 115226 are very close 
to those of known roAp stars. The phase shifts between radial 
velocity variation of different spectral lines suggest the presence 
of running pulsation wave, first propagating in the atmospheric 
layers where hydrogen lines form and then reaching the heights 
where REEs are enhanced. This behaviour is observed in many 
other roAp stars (Ryabchikova et al. 2007b). 

Recently pulsational radial velocity variations at the 
level below 100 ms" 1 were reporte d for HP 13790 9 
dHatzes & Mkrtichianl \2Q04, HD 1 161 14 dElkin et all l2005bl) 
HD 154708 (iKurtz et al.l2006h and HD 75445 (iKochukhov et all 
2007a). In contrast, our discovery of oscillations in HD 1 15226 
adds to the list of roAp stars with much higher amplitudes. In 
fact, given the low probability of obtaining observations at the 
rotation phase of maximum pulsational variability, it is quite 
possible that HD 115226 has one of the highest radial velocity 
amplitudes among all roAp stars. 

Rapid rotation and high pulsat ional ampli t ude m akes 
HD 11522 6 similar to HP 83368 dKochukhovl 12006) and 
HD 99563 dElkin et al.l2005ah . Like these two well-studied roAp 



stars, HD 1 15226 is an obvious candidate for monitoring pulsa- 
tions over complete rotation cycle. We propose to obtain such 
observations both in photometry (to determine rotation period 
and to study the pulsation frequency spectrum) and in spec- 
troscopy (to derive distribution of chemical spots and to study 
the topology of pulsation velocity field). 

The discovery of strong variability in HD 115226 hints 
at the existence of a relationship between the stellar rota- 
tion and excitation of pulsations. The group of roAp stars 
which show radial velocity amplitudes in Ndm at the level 
of 1 kms" 1 and above now includes HD 60435, HD 83368, 
HD 99563, HD 115226, HD 12932 and HD 19918. The first 
four stars all have v e sin/>10 kms _I and P 10t <8 d. Rotation 
periods are unknown for the two latter stars, but cannot be 
long given v e sin i - 3.0-3. 5 kms -1 measured for these objects 
(Ryabchik ova et af]|2007bl) . Thus, we come to the conclusion 
that the high- amplitude roAp stars are markedly different in their 
rotation properties from roAp stars with lower pulsational am- 
plitudes, many of which have rotation periods of the order of 
years and even decades. In the light of this finding, the problem 
of the influence of rotation on excitation of p-modes in roAp 
stars deserves careful theoretical consideration. 
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